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Abstract. We investigate an array of identical phase oscillators non-locally coupled 
without time delay, and find that chimera state with two coherent clusters exists 
which is only reported in delay-coupled systems previously. Moreover, we find that 
the chimera state is not stationary for any finite number of oscillators. The existence 
of the two-cluster chimera state and its time-dependent behaviors for finite number 
of oscillators are confirmed by the theoretical analysis based on the self-consistency 
treatment and the Ott-Antonsen ansatz. 
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1. Introduction 

An array of identical oscillators has been used to model a wide range of systems, such 
as neural networks, convecting fluids, laser arrays and coupled biochemical oscillators. 
These systems exhibit rich collective behaviors including synchrony and spatiotemporal 
chaos [H El [31 H]. Most of the earlier theoretical works on these systems assume 
either local coupling (nearest-neighbor interactions) or global coupling (infinite-range 
interactions); a third type named non-locally coupling began to be explored in the 
past years, which is somewhere between local coupling and global coupling. In non- 
local coupled systems, the oscillators interact with all others and the strength between 
oscillators varies with the distance between them. 

Chimera state is a spatiotemporal pattern in which some of the identical oscillators 
are coherent and synchronous while others remain incoherent [HI EJ El E]. They 
usually appear in systems with non-local coupling and could only be built for proper 
initial conditions [10] . Chimera state does not relate to the partially synchronized states 
observed in populations of nonidentical oscillators with dispersive frequencies in which 
the splitting of the population roots in the inhomogeneity of the oscillator themselves 
and the intrinsically fastest or slowest oscillators remain desynchronized. Its emergence 
cannot be ascribed to a supercritical instability of the spatially uniform oscillation, 
because it occurs even if the uniform state is stable. Chimera may be a paradigm to 
study unihemispheric sleep in neurology which states a fact that many creatures sleep 
with only half of their brain while the other half is still active at the same time. [H] . 

In the year 2002, Chimera state was first reported by Kuramoto and his colleagues 
[3 [6] when simulating the non-locally coupled complex Ginzburg-Landau equation. 
They showed that identical oscillators with non-locally symmetrical coupling could self- 
organize into chimera states. Soon, spiral wave chimera [121 03] was discovered in 
two-dimensional arrays of non-locally coupled oscillators. In succession, Abrams and 
Strogatz (7J found an exact solution for this state in a ring of phase oscillators coupled 
by a cosine kernel. Recently, two interesting findings on chimera state are reported. 
Firstly, in the study of the non-locally coupled oscillators with time delay [HI IT5] . 
clustered chimera state that has spatially distributed phase coherence separated by 
incoherence with adjacent coherent regions in antiphase, was found. The observed 
clustered chimera state in these systems is stationary and its pattern does not change 
with time. Secondly, Abrams and Strogatz [16] found a breathing chimera state in 
a model consisting of two interacting subpopulations of oscillators. Pikovsky and 
Rosenblum [17J considered oscillators ensembles consisting of several subpopulations of 
identical units, with a general heterogeneous coupling between subpopulations, through 
which they acquired quasiperiodic chimera states. Laing [T8l [19] summarized chimera 
states in several heterogeneous networks of coupled phase oscillators, in the mean time, 
he analyzed chimera state applying the Ott-Antonsen ansatz [20] in one-dimensional 
and two-dimensional systems. He pointed out that, in one-dimensional system, when 
parameter heterogeneity is introduced, a breathing chimera state exists. 
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In this work, we study a one-dimensional array of non-locally coupled identical 
phase oscillators. We find that, in the absence of time delay, a two-cluster chimera 
state could exist. We also find that, in the absence of parameter heterogeneity, the two- 
cluster chimera state is not stationary but oscillating. Different from Laing's results |20j, 
we find that the oscillation of the two-cluster chimera state only exists for the system 
with a finite number of oscillators. Both the two-clustered chimera state and finite size 
oscillations of the chimera state in this model are analyzed based on the Ott-Antonsen 
ansatz. 



2. Model 



The array of non-locally coupled phase oscillators can be described in a concise form as 



d(j> 

m =UJ 



/ G(x-x')sm[(/)(x,t) -<p(x\t) +a]dx'. (1) 

J —k 

Here, <p(x, t) is the phase of the oscillator at position x at time t. The space variable 
x is in the range [— k, k] (0 < k < n). The periodic boundary condition is imposed 
for k = 7r, otherwise the no-flux boundary condition is imposed on the system, u is 
the natural frequency (same for all oscillators), which plays no role in the dynamics. 
Without losing generality, we can set u — 0. The angle a (0 < a < |) is a tunable 
parameter. The kernel G(x — x') provides non-local coupling between oscillators. G(x) 
is non- negative, even, decreasing with \x\ along the array, and normalized to have unit 
integral. Following Abrams and Strogatz [TJ, we make use of the cosine kernel 

G ^ = ■ ir (l + Acosx) (2) 

2(k + A sin k) 

where < A < 1. When k = it, Abrams and Strogatz found a chimera state with 
only one coherent cluster [7J. However, we find a novel chimera state which has two 
coherent clusters and is oscillating for any finite number of oscillators. As mentioned 
above, clustered chimera state is only found in the time-delay coupled systems [HI [15] 
and oscillating chimera state is found in the systems with subpopulations [HH |T7] or the 
systems with parameter heterogeneity [i.e., a = a(x)] [T8] . 

Let Q denotes the angular frequency of a rotating frame whose dynamics are 
simplified as much as possible, and let 6 = <p — Qt denotes the phase of an oscillator 
relative to this frame. The key idea behind the analysis of chimera stats is the 
introduction of a mean-field-like quantity, namely, a complex order parameter i?e* e 
[21 El [7J which is defined as 

R{x, t)e l@(x ' t] = I* G(x - x'y e{x,,t) dx' (3) 



—k 



Then Eq (1) becomes 



f)0 

— = u-n-Rsm(6-Q + a). (4) 

For stationary state, R and 9 are time-independent and only depend on space 
variable x. Let A = u> — Q where Q is the angular velocity for the oscillators in coherent 
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Figure 1. (color online) (a) Phase pattern for two-cluster chimera state when the 
steady state is reached. Eq. (1) is integrated using the Runge-Kutta method with fixed 
time step dt = 0.1 with oscillator number N — 256, k — it, ft — 0.10 and A — 0.995. 
(b) The triangle (green) symbol is the modulus R of the complex order parameter, the 
square (black) symbol is the distribution of (0{x)) of individual oscillators averaged 
over 200 time units. The circle (red) symbol is the fluctuation a 2 (x) of 9(x). (c) three 
types of distribution of the phase 6 of the order parameter at different time, (d), 
(e) and (f) show the contour graphs of R(x), Q(x)/n and 8(x)/ir, respectively. The 
horizontal axis is position x and vertical axis is time t. 



regions, the oscillators with A < R are in coherent regions and are phase-locked to 

■R' 



6 = arcsin(f) + 9 - a El U\. 



3. Simulate and results 

For parameters — | — a — 0.10, A = 0.995, k = tc, the system with iV = 256 phase 
oscillators could evolve to a two-cluster chimera state under the initial conditions as 
follows: 

where r is a random variable from a uniform distribution of [—0.5, 0.5]. As shown in 
Fig. 1(a), there exist two coherent clusters in which all oscillators are synchronized. 
Oscillators in the same cluster are nearly in phase yet in different clusters are in 
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Figure 2. (color online) (a) Phase pattern for two-cluster chimera. The parameters 
N = 256, k = 0.8tt, /3 = 0.10 and A = 0.995. (b) The triangle (green) symbol 
is the modulus R of the complex order parameter, the square (black) symbol is the 
distribution of (6(x)) of individual oscillators averaged over 200 time units. The circle 
(red) symbol is the fluctuation a 2 of 9(x). (c) the distribution of the phase O of order 
parameter, (d), (e) and (f) show the contour graphs of R(x), Q(x)/tt and 9(x)/ir, 
respectively. 

antiphase. On the other hand, the oscillators between the two coherent clusters are 
de-synchronized and their phases are randomly distributed in [— tt, 7r]. Then we consider 
three quantities characterizing a chimera state: the modulus R of the complex order 
parameter at an arbitrary time, the angular velocity (0(x)) averaged in a time interval 
of 200 units, and the fluctuation of the instantaneous angular velocity which is defined 
as cr(x) = J ((9(x) — (9(x))) 2 ). The quantities against the locations of oscillators are 
presented in Fig. 1(b). Clearly, there are two plateaus on the curve of (0(x)) which refer 
to the coherent clusters in the chimera state. The zero a in the coherent clusters means 
that the oscillators in the coherent clusters all move on the same instantaneous angular 
velocity. Further, nonzero a outside the coherent clusters refers to the fluctuation 
of angular velocities for the oscillators outside the coherent clusters and indicates 
desynchronization. Fig. 1(b) reveals two features on R for the two-cluster chimera state. 
Firstly, the oscillators can be divided into two domains which join at the minimum of 
R and the curve of R against the locations of oscillators does not show symmetry about 
the minimum of R. Further explorations show that R is a function of time. As shown 
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in Fig. 1(d) where the spatio-temporal evolution of R is featured, R is oscillating in 
each coherent cluster. Especially, when R reaches its minimum in one domain, R in the 
other one reaches its maximum. Secondly, the boundaries of the coherent clusters are not 
determined by the condition of A = R(x) and the coherent regimes are narrower than 
those expected according to A = R(x) in most of the time, which are different from the 
stationary chimera state. The reason for this observation roots in the time-dependent 
order parameter R. Actually, for a forced phase oscillator obeying 9 = uj + R(t)sin6, 
the synchronization of the phase oscillator by the force just requires the phase of the 
oscillator to be confined within (0, 2ir) but not to a fixed value, which leads the onset 
of the synchronization of the oscillator not to obey the condition of A = R(t) and the 
onset of synchronization strongly depends on the details of the functional form of R(t). 
To be noted, even though R is time-dependent, the oscillators in the coherent clusters 
still have the same angular velocity which does not fluctuate as exhibited by zero o. 
The spatiotemporal evolution of R in Fig. 1(d) shows another feature: the two-cluster 
chimera state displays an irregular motion along the ring, for example, the locations of 
coherent clusters vary with time. Similar phenomenon is also observed for the chimera 
state with a single cluster [2TJ. Furthermore, the snapshot and the time evolution of G 
presented in Fig. 1(c) and (e) show that is almost uniform in each domain except for 
those near the junction between domains and there is a phase difference of ir for B in 
different domains. To be stressed, the features on O revealed by Fig. 1(c) and (e) could 
be used as a more general scheme for the initial conditions to generate an oscillating 
two-cluster chimera state. That is, the initial conditions in Eq. [5] could be changed to 
be 0(x, 0) = 7r for x < and (f>(x, 0) = for x > 0. To get a better illustration, we 
present the evolution of 8(x) in Fig. 1(f) which shows that, resulting from the oscillation 
of R, the territories of the coherent clusters alter with time. Furthermore, it should be 
pointed out that, though a two-cluster chimera state could be generated for the system 
characterized by Eq. [1] in the absence of time delay under proper initial conditions, we 
have not detected other chimera states with more than two coherent clusters. 

The two-cluster oscillating chimera state can exist when k ^ tt. In comparison with 
the case with k = ir where oscillators locate on a ring, here the oscillators locate on a 
chain. We take k = 0.8ir as an example. The results are given in Fig. [2j The differences 
with those in Fig. 1 are that the pattern of two-cluster chimera state becomes stationary 
in space and, no matter where the chimera state is initialized, it will adjust its pattern 
to be symmetrical about the center of the chain where the minimum of R appears. 

The analysis above are made on the system with N = 256. One question is how 
the observed two-cluster oscillating chimera state depends on the number of oscillators. 
For this aim, we focus on one domain and record the maximum value of R in this 
domain at any time instance (we denote it as R max )- The time- dependent behavior 
of the two-cluster chimera state can be reflected by the time evolution of R m ax- For 
example, a constant R max indicates a stationary chimera state, otherwise an oscillating 
one. To avoid the influences induced by the irregular motion of the chimera state along 
the ring in the system with periodic boundary condition, we let k = 0.8tt. The results 
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Figure 3. (color online) The maxima of R in one domain varies with time for different 
number of oscillators, (a, b, c, d) are for N = 128, 256, 512 and 1024, respectively. 
Other parameters: A = 0.995, (3 — 0.1 and k = 0.8ir. The (red) lines denote the 
growth of R in the transient and the upper or lower bound of the maxima of R in the 
steady state. The inset in (d) shows the amplitude of the oscillation of R m ax, AR max , 
against the system size N, which indicates a power law of AR max ~ N~ - 5 . 



for different N are presented in Fig. |3j One remarkable feature revealed by the figure 
is that, before the oscillating chimera state is established, the system first evolves to a 
two-cluster chimera state which looks like a stationary one due to the weak oscillation 
of Rmax- However, the "stationary" state is not stable and its instability leads to the 
appearance of an oscillating chimera state. Interestingly, the stationary value of R ma x 
is independent of the size of the system and the time consumed for the system to build 
an oscillating state becomes much longer as the number of the oscillator N increases. 
Another feature in FigJ3 is that larger N seems to weaken the oscillation of R max in the 
two-cluster oscillating chimera state, which is prominent by the comparison between Fig. 
El(a) (N = 128) and (b) (N = 256). Since the transient time to build an oscillating two- 
cluster chimera state for large N becomes extremely long, we just give a rough estimate 
based on the data presented in the inset of Fig. EJ^d) that the oscillation amplitude 
AR max ~ iV~ a5 , which means that the two-cluster oscillating chimera state becomes a 
stationary one in the thermodynamic limit, which is different from the Laing's results 

ua. 
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Figure 4. (color online) The modulus R and the phase of the order parameter and 
A by solving of Eq. ^ via an iterative scheme as presented in the context, (a) The 
square (black) symbol denotes the distribution of R and the circle (red) the value of 
A. (b) Distribution of O. Other parameters: A = 0.995, /3 = 0.1 and k = it. 



4. Analysis 

The above results can be understood theoretically. First, the two-cluster stationary 
chimera state in the thermodynamic limit can be explained in terms of Kuramoto- 
Battogtokh self-consistency equation [5] as follows: 

•« r - a , n A - \/a 2 - R(x') 2 

R( x y e(x) = e t(3 / G(x - x')e ie{x ) x V — dx' (6) 

J-TT R(X') 

note that there are three unknown quantities (the real- valued functions R(x), Q(x) and 
the real number A) in terms of the assumed choices of /3 and the kernel G(x — x'). We 
take k = n as an example. Considering the patterns in Figs. 1(b) and (c), the modulus 
R and the phase of the order parameter approximately satisfy 

R(x + tt) = R{x) 

Q(x + n) = Q{x)+7i. 1 ' 

We substitute Eqs. © and (f5]) into Eq. (JSJ and get 

R(x)e ie{x) = e i/3 [° — [1 + Acos(x - x')e ie ^H(x')dx' 

+ ^ r — t 1 + A cos ( x - x')}e ie{x ' ] H(x')dx', (8) 
jo 2% 

where H(x) = (A — ^A 2 — R(x) 2 )/R(x). Under the transformation x' — n — > x' in the 
second term of the right hand side of Eq. flS]), the self-consistency equation changes into 

R( x )e ie ^ = e i(3 [° AC0S{x ~ X,) e ie ^H(x')dx'. (9) 

J-TT 7T 
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Figure 5. (color online) Phase pattern for two-cluster chimera when the steady state 
is reached. Eq. p7|) and (fT5|) is integrated using the Runge-Kutta method with fixed 
time step dt = 0.005 and the oscillator number N — 256. (a, c) are the modulus R and 
the phase Q/tt of the complex order parameter with k — tt. (b, d) are the modulus R 
and the phase Q/tt of the complex order parameter with k = 0.87T. Other parameters: 
13 = 0.1, A = 0.995. 



To solve Eq. (jHJ), we first determinate the value of A. Because Eq. fl2]) is left unchanged 
by any rigid rotation Q(x) — > Q(x) + Bo, we can specify the value of Q(x) at any point 
x we like. We set O(f) = 0. Now we can get A. Then we take R(x) and Q(x) 
obtained from the dynamical simulations as initial guesses and use an iterative scheme 
to determinate R(x) and Q(x) in function space, behind which the idea is that the 
current estimates of R(x) and Q(x) can be entered into the right-hand side of and 
used to generate the new estimates appearing on the left-hand side. Figure @] shows 
the results obtained from Eqs.(9) and (2). To be stressed, without the requirement of 
Eq.(7), the self-consistency equation for any finite system always yields to a one-cluster 
chimera state, which also evidences that the stationary two-cluster chimera state here 
is not stable. 

From Fig. HI we could notice the order parameter of stationary state, which has 
a little difference from Fig. [Tfb) and (c) in the vicinity of the junction between two 
domains. The transition of the former case performances like a step function while the 
latter one is continuous. The difference probably originates from the finite size effect in 

Fig. m 

The oscillating characteristic could be interpreted with the assistance of the Ott- 
Antonsen ansatz [TH [20]. Following the line in [161 El EEE], we assume that there is 
a probability density function f(x,co,9,t) characterizing the state of the system. This 



10 



0.4 



P 0.2 



0.0 



-■-0.995 
-•-0.95 




. -f. \ 






— / \ 

{ 




■ 

/ 

■ 

/ 


A 


■ 

/ 

■ 

/ 

■ 


/ 

• 

/ 

• 






—¥- « 



0.6 



0.8 

k 



1.0 



Figure 6. (color online) The boundary of the two-cluster chimera state plotted on 
k — P parameter plane. Below the curves, the two-cluster chimera state is stable. The 
square (black) symbol and circle (red) symbol are plotted at the parameter values 
A = 0.995 and A = 0.95, respectively. 



function satisfies the continuity equation |T6l [T7l HH [19] 



where 



v = iu — G(x — x') 
-k 



oo PIT 



(10) 



sin(fl - 6' + a)f(x', oo, 6', t)d6'dujdx' (11) 



OO «/ —TV 



with 9 = 9(x) and & = 9(x f ). The complex order parameter can be formulated as 

/k roo rn 

G{x-x') / / e id f(x',u,9',t)d9'dujdx'. 

In terms of complex order parameter Z, Eq. (Til]) can be rewritten as 

v = u-^[Ze- i <- e -® + Z*e^} 

where Z* denotes the complex conjugate of Z and f3 = | — a. Following Ott and 
Antonsen [181 EQl , we have 



(12) 



(13) 



f(x,u,9,t) 



2tt 



{i + J2Mx,u,ty e r + c.c.}} 



(14) 



n=l 



where c.c. is the complex conjugate of the previous term and g(u)) is the distribution 
of natural frequency. In this work, we assign all oscillators a same natural frequency 
{u = 0), so g(u) = 6(u). Substituting Eqs.JISJ and (dU) into Eqs. and (JT2J, we 
obtain 

da(x,u,t) i 



dt 



■\Z*e- ip + Ze if3 a 2 



(15) 
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Figure 7. (color online) Two-cluster chimera under the condition of different types of 
nonlocal coupling kernel G(x). The top panels are the phase patterns for two-cluster 
chimera and the bottom are the modulus R (black) and the phase (red) of the order 
parameters, (a, b) is computed under the assumption that G(x) takes exponential 
form with parameters: A = 4, k = 1 and j3 = 0.10 while G(x) in (c, d) is steplike 
function with parameters: d — 0.5, k = 1 and (3 — 0.10. All above take no-flux 
boundary condition and the initial conditions are as following: 4>{x, 0) = for x < 
and 4>{x, 0) = 7r for x > 



Z=j k G{x-x')l g(co)a*(x,u,t)dujdx'. (16) 



-k 

Letting a(x,t) = a(x,0,t), we have 



^A = h Z * e ^ + Ze^a 2 } (17) 

k 

G(x-x')a*(x',t)dx'. (18) 

-k 



By numerically simulating these two equations, we have the time evolutions of R(x) 
and Q(x). The results for k = n and k = 0.8tt are presented in Fig. [5j respectively. 
Clearly, the main features in Figs. (1) and (2), such as the oscillating nature of the 
two-cluster chimera state, the movement of the pattern (or the frozen pattern) of R and 
for k = 7T (or for k = 0.871"), and the uniform distribution of B(x) in different domains, 
are reproduced in Fig. [51 
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Using Eqs. (I17p and (1181) . we may probe into the regime for the existence of the 
two-cluster chimera state on k — (3 parameter plane. The results are presented in Fig. [6] 
at A = 0.995 and A = 0.95. As shown in this plot, the two-cluster chimera state is not 
favorable at large /?, and either small k or large k tends to be harmful for the two-cluster 
chimera state. And the smaller A is, the smaller the domain of two-clustered chimera 
exists, if A < 0.8, the domain does not exist any longer. Beyond the stable regime for 
the two-cluster chimera, i.e., above the curves in Fig. [61 the initial two-cluster chimera 
state tends to become a normal chimera state with only one coherent cluster. The two- 
cluster oscillating chimera state may be detected in coupled oscillators with other types 
of non- locally coupling kernel G(x) which could be exemplified by two cases. In the first 

A —A\x\ 

one, G{x) takes an exponentially decaying function, G exp {x) = 2 (i-e- Ak ) • ^ n ^ ne °^er 
case, G(x) takes a steplike function, G step (x,d) = ^ for |x| < d and G step (x,d) = for 
| a; | > d. Figure 7 shows the results for the systems with these two types of non- locally 
coupling and, clearly, the two-cluster oscillating chimera states are reproduced. 

5. Conclusion 

In summary, we study a one- dimensional system consisting of non-locally coupled phase 
oscillators, which is a prototype for studying chimera states. By numerically simulating 
this simplest system, we find the existence of a two-cluster oscillating chimera state in 
the absence of time delay coupling and parameter heterogeneity. The numerical results 
are confirmed by the theoretical analysis based on the self-consistency treatment and 
the Ott-Antonsen ansatz. 
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